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We consider the cross section in Fourier space, conjugate to the outgoing hadron's transverse mo- 
mentum, where convolutions of transverse momentum dependent parton distribution functions 
and fragmentation functions become simple products. Individual asymmetric terms in the cross 
section can be projected out by means of a generalized set of weights involving Bessel functions. 
Advantages of employing these Bessel weights are that they suppress (divergent) contributions 
from high transverse momentum and that soft factors cancel in (Bessel-) weighted asymmetries. 
Also, the resulting compact expressions immediately connect to previous work on evolution equa- 
tions for transverse momentum dependent parton distribution and fragmentation functions and to 
quantities accessible in lattice QCD. Bessel-weighted asymmetries are thus model independent 
observables that augment the description and our understanding of correlations of spin and mo- 
mentum in nucleon structure. 
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1. Weighted asymmetries 



In the factorized picture of semi-inclusive processes, where the transverse momentum of the 
detected hadron P/,j_ is small compared to the photon virtuality Q 2 , transverse momentum depen- 
dent (TMD) parton distribution functions (PDFs) characterize the spin and momentum structure 
of the proton [jl], [| f| f§ Q g, [7p. At leading twist there are 8 TMD PDFs. They can be studied 
experimentally by analyzing angular modulations in the differential cross section, so called spin 
and azimuthal asymmetries. These modulations are a function of the azimuthal angles of the final 
state hadron momentum about the virtual photon direction, as well as that of the target polarization 
(see Fig. |] and e.g., Ref. [||] for a review). TMD PDFs enter the SIDIS cross section in mo- 
mentum space convoluted with transverse momentum dependent fragmentation functions (TMD 
FFs). However, after a two-dimensional Fourier transform of the cross section with respect to 
the transverse hadron momentum P/jj_, these convolutions become simple products of functions in 
Fourier b^-space. The usefulness of such Fourier-Bessel transforms in studying the factorization 
as well as the scale dependence of transverse momentum dependent cross section has been known 
for some time [Q, 10, 11, 12, 13]. Here we exhibit the structure of the cross section in b^-space 



and demonstrate how this representation results in model independent observables [ 14 ] which are 



generalizations of the conventional weighted asymmetries fll5|, |6|, Additionally, we explore the 
impact that these observables have in studying the scale dependence of the SIDIS cross section at 
small to moderate transverse momentum where the TMD framework is designed to give a good 
description of the cross section. In particular we study how the so called soft factor cancels from 
these observables. The soft factor g|, |l7|, [IT], @ HJ is an essential element of the cross section 
that arises in TMD factorization expressions [^, [R], |Tl], [13| |. It accounts for the collective effect 
of soft momentum gluons not associated with either the distribution or fragmentation part of the 



process and it is shown to be universal in hard processes [17]. 

The concept of transverse momentum weighted Single Spin Asymmetries (SSA) was proposed 
some time ago in Ref. [|[ [7]]. Using the technique of weighting enables one to disentangle in a 
model independent way the cross sections and asymmetries in terms of the transverse momentum 
moments of TMDs. A comprehensive list of such weights was derived in Ref. [0] for semi-inclusive 
deep inelastic scattering (SIDIS). In SIDIS and Drell-Yan scattering, proofs of TMD factorization 
contain an additional factor, the soft factor [ITS, 13, 17, 11, 12, 18]. At tree level (zeroth order in 



as) the soft factor is unity, which explains its absence in the factorization formalism considered 
for example in Ref. [^J. Consequently it is also absent in tree level phenomenological analyses 
of the experimental data (see for example Refs. [|l^, 21, 22]). In principle, the results of tree 



level analyses at different energies cannot compared. For a correct description of the energy scale 
dependence of the cross sections and asymmetries involving TMDs, the soft factor is essential to 



include nTTl 12]. However, it is possible to consider observables where the soft factor is indeed 



absent or cancels out. These are precisely the weighted asymmetries [14]. 

We focus on the Sivers asymmetry |2p. With a general |P/,jJ -weight, this asymmetry can be 
written as [0] 



2f d\P h x\d(j) h dQ s wi(\Ph±\) sin(^-^s) (dc(<$> h ,<$> s ) -do^fc + n)) 



Jd\P h ±\d(j) h d(j) S wo(\P h ±\) (da(^ h ,^ s ) + da(^ h ^ s + 7t)) 



(1.1) 
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In the numerator, the angular weight sin(0/, — 0j) projects out the structure function F^p' ^ [||] 
from the cross section, while the weight-wi (|P/u_|) leads to a "deconvolution" of the structure 
function into a product of the first p^-moment of the Sivers function /j^ and the lowest moment 
of the unpolarized fragmentation function Dj [^]. In a similar manner using the weight wo(|P^_l|) 
the unpolarized structure function Fuu.t is written zeroth p^-moments, f\ (x) and Df*\ 



In []14p, we show that in TMD factorization the soft factor cancels in the asymmetry due to the 
"deconvolution" achieved by appropriate |P/u_| -weighting. We demonstrate that it is natural to em- 
ploy Bessel functions as weights, w n J n (\Ph±\^T), where SSj is the Fourier conjugate variable 
of |Paj_|. This generalized weighting procedure addresses a problem related to the perturbative tail 
of TMDs when we weight with conventional weights [7p, w n |Pftj_|' ! . Using Bessel functions 
as weights, the respective integrals become convergent while the "deconvolution" property and soft 
factor cancellation are preserved. However, it is important to stress that while the integrals are con- 
vergent, the scale dependence of the resulting functions, and consequently also the Q 2 dependence 
of the Bessel-weighted asymmetries remains to be studied. 

1.1 Bessel-weighted asymmetries 

To deconvolute or convert the convolutions of TMD PDFs and TMD FFs in the SIDIS cross 
section into products, one can perform a multipole expansion and a subsequent Fourier transform 
of the cross section with respect to the transverse components P/,j_ of the hadron momentum. In 
general, one can write a transverse momentum dependent cross section a(\Ph±\^h) as a two- 
dimensional multipole expansion of the cross section in Fourier space [14]. The n th harmonic in tpt, 



is given by the n th Bessel function of the first kind /„. With such definitions, the relevant terms of 
the SIDIS cross section can be written as 

da a 2 y 2 f^ + Y 2 \ f" d\b T \ 



dxdyd(j>sdzd(j)i l \Ph±\ d\Ph±\ xyQ 2 {l — e)\ 2x ) Jo (In) 

x j/ () (|b r ||P/,J)iW + |S±|sin(^-^)/i(|br||P*j.|)F^-* l) + ... j, (1.2) 

where the ellipsis represent 16 more terms. Note that there is only a finite number of multipoles 
in the SIDIS cross section; the Bessel function of highest order is Jt,. Here we show only two 
terms in the cross section and omit for now regularization parameters needed beyond tree level, see 
Ref. [ |l4| ] and references therein for more details. Introducing the Fourier-transformed TMDs and 
fragmentation functions, 

f(x,b 2 ) = J d 2 p T J** f(x,p 2 ) =2n £ d\p T \\p T \J (\b T \\p T \)f(xy T ) , 

and the derivative (or bj moment) 

2 27T f°° \d 2 I 

f^(x,b 2 T )= - w d b 2f(x,b 2 ) = w j o rf|pr|gU(|br||pr|)/(*,p 2 r), (1.3) 

the structure functions in Fourier space in Eq. ( |L"2| ) are given by 

Fuu.t = x^el ~ff )a {x,z 2 b 2 ) D? )a (z,b 2 ) S(b 2 ) H UVJ {Q 2 ) , (1.4) 

a 

F$t* s) = -x^MzM fti l) \x,z 2 b 2 )bf\z,b 2 ) S(b 2 )H UT . T (Q 2 ) , (1.5) 
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Figure 1; Kinematics of the SIDIS process. The in- and out-going lepton momenta are / and /', respectively. The momentum transfer 
is q. The target nucleon carries momentum P and its transverse spin components are labelled S± . The momentum of the measured 
hadron P/, has transverse components P/,x, which define an angle 0/, with the lepton plane. 



where /j and /{y are the Fourier transformed unpolarized and Sivers TMD PDFs respectively, 
and D?> is a Fourier transformed fragmentation function. We have used the kinematic variables 
Q 2 = -q 2 , M 2 = P 2 , x w x B = Q 2 /P ■ q, y = P ■ q/P ■ I, and z ~ z h = P • Ph/P ■ q and assume 
|P/,_l| <C z.Q- The sum £ a runs over quark flavors a and e a is the corresponding electric charge of 
the quark. In contrast to the tree-level equation [Jsj, |i4[ |, we include here an explicit soft factor S(bj) 
(in Fourier space) and a scale dependent hard part H{Q 2 ), as in reference [11, 12]. For brevity we 
suppress the dependencies on a renormalization scale /J. and on rapidity cutoff parameters (e.g., 
£, f , p in [11, 12, 0]). It is now clear that the cross section Eq. (L2) is a multipole series, 



and that projection on Fourier modes in polar coordinates fa, |P/u_| will give access to the right 
hand sides of Eqs. ( |1.4f ) and (|T3|). Calculating the weighted asymmetry Eq. (1.1) with weights 
Wi(|Paj_|) = 2J l (\F h± \.^ T )/zM^T and w (|P*j.|) =M\Phx\^T) thus yields 

A ^§^Mh-^) ( . _ , _ -2I u , u 2 // c/r , r (G 2 )/ 1 L r (1)fl (x,z 2 ^ ) ^2 )5 (0) fl(z ^2 ) 



(1.6) 

Due to the "deconvolution" of the structure functions in the weighted asymmetries, and universality 
of the soft factor, S cancels in the numerator and the denominator. 

Further, note that ,%j enters the weights wq and wi as a free parameter that we can scan over 
a whole range in order to compare the transverse momentum dependence of the distributions in the 
numerator and denominator relative to each other (in Fourier space). At the operator level, £%t{= 
|br|) controls the space-like transverse distance between quark fields in the correlation functions 
we measure. The Bessel-weighted asymmetries are a natural extension of conventional weighted 
asymmetries [^, [7|] with weights w\ proportional to powers of |Pju_|. Indeed, in the limit S$j — > 0, 



Eq. (1.6) results in the often encountered special case for the SIDIS cross section 



ip 



Ajf Sm ^ vs '(x,z,y) = -2 ^ a " U1 ' iy ~ > Jl ' w , (1.7) 

L a elH UUtT (Q 2 )fi° )a (x)D? )a (z) 

where formally f[°\x,0) = f[°\x), fi T W (x,0) = f^ l \x), and fij 0) (x,0) = Df\x) are 

/»(x,0) =f^(x) = f d 2 PT (jjjL^ f(x,p 2 ). (1.8) 

However we caution that these moments are not well-defined without some regularization. It is 
therefore safer to study Bessel-weighted asymmetries at finite SS T , where the Bessel functions 



^_ ,L a e 2 a H UT , T (Q 2 )ft T [1)a (x)Df )a (z) 
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suppress contributions from large transverse momenta. Without further regularization, the inte- 
grals defining the moments /j^! -(x) and f^(x) are ill-defined due to the asymptotic behavior 
/ir(*>Pr) °* 1/Pr. fii^Pr) « VPr at lar § e Pr (see |§|l). 

1.2 Cancellation of soft factor at the level of matrix elements 

In a similar manner to the discussion above, we now consider the soft factor cancellation in the 
average transverse momentum shift of unpolarized quarks in a transversely polarized nucleon for a 
given longitudinal momentum fraction x [|4|] defined by a ratio of the pj -weighted correlator [||] 



TU 



Jd 2 p T p y ^r + \x,p T ,P,S,^^,p) 



M 



flT 'fan ,C,P) 



5±=0,S r =(l,0) f\ O^M^CiP 



(1.9) 



where /jy and /j are the moments defined in Eq. (1.8). Obviously, the average transverse mo- 
mentum shift is very similar in structure to the weighted asymmetry Eq. ( |l.6| ). While the weighted 
asymmetries are accessible directly from the P h ± -weighted cross section, the average transverse 
momentum shifts are obtained from the pj -weighted correlator and could in principle be accessi- 
ble from weighted jet SIDIS asymmetries. Therefore we generalize the above quantity, weighting 
with Bessel functions of |pr|. In particular, we replace 

2Ji(\p T \r 



sin(0 p - </) 5 ) 



(1.10) 



where 05 
thus 



p y = \p T \ sin(0 p ) — 
for the choice Sj = (1,0) in Eq. (l.S). The Bessel- weighted analog of Eq. (1.9) is 



{Py{x)) T U = 



Jd\p T \\p T \Jd<p / Ml %f T h m($ p -fc) <i>^(x,p T ,P,S,li 2 ,C,p) 



fd\p T \\p T \Jd(l> p Jo(\PT\&T) 
f^\x,3S 2 T ;V- 2 X,P) 



&+)iyi(x,p T ,P,S,Li 2 ,Z,p) 



M 



(1.11) 



fi°\x,m 2 ;n\C,p) ' 

where the correlator <!>(+) is given in p5|]. Again, the soft factors cancel where the independence 
of the soft factor on v-b/Vv 2 is crucial [14]. Further, weighting with Bessel functions at various 



lengths &t thus allows us to map out, e.g., ratios of Fourier-transformed TMDs p4|]. In the limit 
3Bj — > 0, we recover Eq. (|L9|), {py(x)}j- v = (p y (x))Tu > which we have thus shown to be formally 
free of any soft factor contribution. Again we caution that the expressions at S3j = are ill-defined 
without an additional regularization. Finally we note that the study of the ratio Eq. ( |1 . 1 1\ > in lattice 
QCD gives direct numerical evidence that non-zero T-odd TMDs are a consequence of the first 



principles of QCD [25] 



1.3 Conclusions 

Rewriting the SIDIS cross-section in coordinate space displays the important feature that struc- 
ture functions become simple products of Fourier transformed TMD PDFs and FFs, or derivatives 
thereof. The angular structure of the cross section naturally suggests weighting with Bessel func- 
tions in order to project out these Fourier-Bessel transformed distributions, which serve as well- 
defined replacements of the transverse moments entering conventional weighted asymmetries. In 
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addition, Bessel- weighted asymmetries provide a unique opportunity to study nucleon structure in 
a model independent way due to the absence of the soft factor which cancels from these observ- 
ables. This cancellation is based on the fact that the soft factor is flavor blind in hard processes, 



and it depends only on b 2 -,^ 2 and the rapidity cutoff parameter p [11, 12 j. Moreover, evolution 
equations for the distributions are typically calculated in terms of the (derivatives of) Fourier trans- 
formed TMD PDFs and FFs. As a result, the study of the scale dependence of Bessel-weighted 



asymmetries should prove more straightforward [126, 27, Eq]. Thus, we propose Bessel-weighted 



asymmetries as clean observables to study the scale dependence of TMD PDFs and FFs. 
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